Thisfull  text  paper  was  peer  reviewed  at  the  direction  of  IEEE  Communications  Society  subject  matter  experts  for  publication  in  theWCNC  2010  proceedings. 


Many-to-Many  Communications  via  Space-Time 

Network  Coding 

Mohammed  W.  Baidas  1 ,  Hung-Quoc  Lai  2  and  K.  J.  Ray  Liu  1 

1  Department  of  Electrical  and  Computer  Engineering,  University  of  Maryland,  College  Park,  MD  20742,  USA 
2  US  Army  RDECOM  CERDEC,  RDER-STW-NA,  Bldg.  245,  Aberdeen  Proving  Ground,  MD  21005,  USA 
(email:  baidas@ieee.org,  hungquoc.lai@us.army.mil  and  kjrliu@umd.edu) 


Abstract — In  this  paper,  the  mutual  cooperative  communication 
between  multiple  nodes  in  a  wireless  network  is  efficiently  achieved 
through  a  novel  concept  of  Space-Time  Network  Coding  (STNC).  Unlike 
the  conventional  point-to-point  cooperative  communications  between 
two  nodes  with  TV  relay  nodes  deployed  in  between,  simultaneous 
transmissions  from  the  different  TV  nodes  acting  as  source/relay  nodes 
are  performed  within  2 TV  time-slots.  In  particular,  the  communication  is 
split  into  two  phases:  1)  Broadcasting  Phase  and  2)  Cooperation  Phase. 
In  the  Broadcasting  Phase,  each  node  broadcasts  its  data  symbol  to  the 
other  nodes  in  the  network  in  its  own  time-slot,  alternatively;  while  in 
the  Cooperation  Phase,  in  each  time-slot,  a  set  of  (TV  —  1)  nodes  transmit 
while  a  single  destination  node  receives  the  other  nodes’  transmissions. 
Specifically,  each  node  employing  a  selective  Decode-and-Forward  (DF) 
cooperative  protocol,  performs  a  linear  combination  of  the  other  nodes 
data  symbols  and  all  the  (TV  —  1)  nodes  simultaneously  transmit  their 
signals  to  a  single  receiving  node;  which  then  performs  joint  multiuser 
detection  to  separate  the  different  nodes’  symbols.  Exact  symbol-error- 
rate  (SER)  expressions  for  arbitrary  order  M-ary  Phase  Shift  Keying 
(M-PSK)  modulation  are  derived.  In  addition,  an  asymptotic  SER 
approximation  is  also  provided  which  is  shown  to  be  tight  at  high 
signal-to-noise  ratio  (SNR).  Finally,  the  analytical  results  confirm  that 
for  a  network  of  TV  nodes,  a  full  diversity  order  of  (TV  —  1)  per  node  is 
achieved  by  the  proposed  STNC  cooperative  communication  scheme. 

Index  Terms — Cooperation,  decode-and-forward,  diversity,  many-to- 
many,  symbol-error-rate,  wireless  network  coding 

I.  Introduction 

NETWORK  CODING  has  recently  emerged  as  an  effective  ap¬ 
proach  for  efficiently  distributing  data  and  increasing  network 
throughput  [1]  [2].  Cooperative  communications  have  also  attracted 
much  attention  in  the  research  literature  due  to  their  achievable 
spatial  diversity  gains,  enhanced  coverage  and  improved  transmission 
reliability  [3]  [4].  In  conventional  relay  networks,  a  set  of  TV  relay 
nodes  are  deployed  between  the  source  and  destination  nodes  and 
the  available  network  bandwidth  is  split  into  TV  +  1  orthogonal 
channels  using  TDMA.  However,  with  the  increase  in  the  number 
of  relay  nodes,  the  conventional  multinode  relay  networks  become 
excessively  bandwidth  inefficient.  Moreover,  the  classical  cooperative 
communication  protocols  are  not  well  suited  for  distributing  infor¬ 
mation  from  one  or  more  source  nodes  to  possibly  many  destination 
nodes,  simultaneously.  This  in  turn  suggests  the  use  of  the  concept 
of  wireless  network  coding  for  exchanging  data  symbols  among 
multiple  cooperative  nodes  over  wireless  networks. 

Recently,  there  have  been  several  research  works  aiming  at 
employing  wireless  network  coding  in  cooperative  networks.  For 
instance,  in  [5],  location- aware  cooperative  wireless  network  coding 
through  the  novel  concept  of  Wireless  Network  Cocast  (WNC)  was 
proposed.  In  particular,  it  was  illustrated  that  with  WNC,  a  reduction 
in  aggregate  transmission  power  and  delay  can  be  achieved  along 
with  incremental  diversity  for  different  relaying  schemes.  In  [6], 
an  algebraic  superposition  of  channel  codes  over  a  finite  field  is 
proposed  to  allow  two  nodes  to  cooperate  in  transmitting  information 


to  a  single  destination.  Bi-directional  relaying  between  two  source 
nodes  through  a  single  relay  node  employing  wireless  network 
coding  has  been  introduced  in  [7].  An  outage  analysis  of  network 
coded  communication  of  multiple  users  with  a  single  destination 
node  through  a  set  of  dedicated  relay  nodes  was  analyzed  in  [8]. 

In  this  paper,  the  merits  of  network  coding  and  cooperative 
diversity  are  exploited  to  allow  TV  nodes  ’’users”  to  exchange  data 
between  each  other  with  the  novel  concept  of  Space-Time  Network 
Coding  (STNC).  In  this  work,  the  Decode-and-Forward  (DF)  relaying 
protocol  [3]  is  studied  within  the  concept  of  STNC.  In  particular, 
the  STNC  scheme  is  based  on  linear  wireless  network  coding  over 
DF  nodes  and  the  communication  is  split  into  two  phases:  1) 
Broadcasting  Phase  and  2)  Cooperation  Phase.  In  the  former  phase, 
each  node  broadcasts  its  data  to  the  other  nodes  in  its  dedicated  time 
slot;  while  in  the  latter  phase,  in  each  time-slot,  a  set  of  (TV— 1)  nodes 
transmit  linearly  network  coded  signals  while  a  single  destination 
node  receives  the  other  nodes’  transmissions.  A  simple  multiuser 
detection  [9]  is  then  applied  at  each  node  to  separate  the  different 
data  symbols  received  from  the  different  nodes. 

Exact  analytical  derivations  of  the  symbol  error  rate  (SER)  perfor¬ 
mance  of  the  STNC  scheme  and  comparative  simulation  results  are 
provided  in  this  paper.  Moreover,  tight  asymptotic  approximation 
at  high  signal-to-noise  ratio  is  also  analyzed  and  the  cooperative 
diversity  order  achievable  with  the  STNC  scheme  is  verified.  It  is 
illustrated  that  with  the  STNC  scheme,  N  information  symbols  of 
all  the  N  nodes  can  be  exchanged  over  a  total  of  2TV  time-slots  (i.e. 
1/2  symbol  per  node  per  channel  use)  as  well  as  achieving  a  full 
diversity  order  of  (TV  —  1)  per  node. 

In  the  remainder  of  this  paper,  the  system  model  and  communi¬ 
cation  phases  are  presented  in  Section  II.  The  multi-source  signal 
detection  and  the  exact  SER  analysis  are  presented  in  Sections 
III  and  IV,  respectively.  In  Section  V,  the  asymptotic  upper-bound 
SER  expression  is  derived  while  the  SER  performance  evaluation  is 
presented  in  Section  VI.  The  conclusions  are  drawn  in  Section  VII. 

II.  System  Model 

Consider  a  wireless  network  consisting  of  TV  nodes  (TV  >  3) 
denoted  as  S2,  •  •  • ,  Sn-  In  this  model,  each  node  is  equipped 
with  only  one  antenna  and  can  act  as  a  source,  relay  or  destination. 
Without  loss  of  generality,  the  Decode-and-Forward  (DF)  coopera¬ 
tion  protocol  is  considered.  The  TV  nodes  are  assumed  to  have  their 
own  information  symbols  as  xi,  X2, . . . ,  xn,  respectively,  and  each 
node  wishes  to  exchange  its  own  data  symbol  with  the  other  nodes. 

The  communication  between  all  the  source  nodes  is  split  into  two 
main  phases,  namely  the  Broadcasting  Phase  and  the  Cooperation 
Phase,  over  a  total  of  2TV  time-slots,  TV  time-slots  each.  During  the 
broadcasting  phase,  source  node  Sj  is  assigned  a  time-slot  Tj  in 
which  it  broadcasts  its  own  data  symbol  to  the  other  nodes  Si  for 
i  e  {1,2,...,  TV}  for  i  ^  j.  That  is,  the  broadcasting  phase  is  an 
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Fig.  1.  Space-Time  Network  Coding  -  Broadcasting  and  Cooperation  Phases  -  N  =  4 
Nodes 

information  exchange  phase,  which  upon  completion,  each  node  Si 
will  have  received  a  set  of  (TV  —  1)  symbols  xi, . . . ,  Xj, . . . ,  xjy  for 
j  7^  i  from  the  other  nodes.  With  respect  to  the  cooperation  phase, 
each  node  acts  as  a  relay  for  the  other  nodes  with  one  node  being  the 
destination  in  each  time-slot.  In  particular,  each  node,  except  a  single 
receiving  node,  forms  a  linearly-coded  signal  from  the  overheard 
symbols  and  transmits  it  to  the  receiving  node,  which  upon  receiving 
the  (TV  —  1)  linearly- coded  signals,  performs  a  multi-source  signal 
separation  to  extract  the  desired  symbol  from  each  node.  To  allow  for 
joint  detection/source  separation  for  each  of  the  linearly-coded  trans¬ 
mitted  symbols  of  the  different  nodes  at  each  receiving  node,  each 
symbol  Xj  is  spread  using  a  signature  waveform  Sj(t).  The  cross¬ 
correlation  between  waveforms  Sj  (£)  and  Si  (t)  over  a  symbol  dura- 
tion  Ts  is  given  by  pjti  =  {sj(t),Si(t))  =  (1  /Ts)  /QTs  Sj(t)s*(t)dt, 
where  it  is  assumed  that  pjj  m  ||sj(t)||2  =  1.  It  is  further  assumed 
that  each  node  knows  the  signature  waveform  of  the  other  nodes 
in  the  network  which  is  required  for  multi- source  detection.  The 
cooperative  STNC  communication  scheme  is  illustrated  in  Fig.  1. 

A.  Broadcasting  Phase 

In  this  subsection,  the  signal  model  for  an  arbitrary  symbol  Xj, 
j  E  {1,2,...  ,2V},  transmitted  during  the  broadcasting  phase  from 
node  Sj  to  the  other  nodes  is  presented.  Node  Sj  broadcasts  it  own 
data  symbol  Xj  to  the  other  nodes  in  its  dedicated  jth  time-slot. 
Thus,  the  signal  received  at  each  node  Si  for  i  ^  j  is  given  by 

PSj  hjiXj Sj(t)  T  rij^(t),  (1) 

where  PSj  is  the  transmitted  power  by  node  Sj  (the  power  allocation 
will  be  discussed  in  a  later  section),  Sj  ( t )  is  the  signature  waveform 
of  node  Sj,  and  rij^{t)  is  the  additive  white  Gaussian  noise  sample 
at  node  Si  due  to  the  signal  transmitted  by  node  Sj.  In  (1),  nj^{t) 
is  modeled  as  a  zero-mean  complex  Gaussian  random  variable  with 
variance  TV0.  Moreover,  hj^  is  the  flat  fading  channel  coefficient 
between  nodes  Sj  and  Si  that  is  distributed  as  a  zero-mean  complex 
Gaussian  random  variable  hj ^  ~  CTV(0,  cr2^),  where  a2^  is  the 
channel  gain.  In  addition,  hj ^  can  be  expressed  as  hj 
where  \hj^\  —  ctj ^  is  the  Rayleigh  distributed  magnitude  as 

.  2a  (  a2  \ 

=  ~2~  exP  I  -~2~  I  j  «  >  (2) 

aj,i  \  <7J,i  J 

with  (j)j^  being  the  phase  response,  uniformly  distributed  over  the 
interval  [— 7r,7r).  Moreover,  it  is  assumed  that  the  receiving  node  Si 
can  perfectly  estimate  the  channel  coefficient  hj ^  from  the  received 
signal  yj^i.  Also,  the  channels  are  assumed  to  be  reciprocal  (i.e. 
hij  =  hj:i)  as  in  Time  Division  Duplexing  (TDD)  systems. 


The  broadcasting  phase  can  be  put  in  matrix  form  as  follows 


Si 

...  Sj 

Sn 

Ti 

~VPs  1X± 

0 

0 

Tj 

0 

y/PsjXj 

0 

Tn 

_  0 

0 

•  s/Ps  N  XN- 

The  detection  of  the  data  symbol  Xj  at  the  node  Si  can  be 
achieved  by  cross-correlating  the  received  signal  Vj,i(t)  in  (1)  with 
the  signature  waveform  Sj{t)  as 


vj,i  =  <%,»(*),  «j(*))  = 


(4) 


where  rij ^  ~  CJ\f(0,  No).  Upon  the  completion  of  broadcasting 
phase  (i.e.  after  TV  time-slots),  each  node  Si  will  have  received  a 
set  of  (TV  —  1)  symbols  {yj,i}f=1  ^  from  all  the  other  nodes  in  the 
network.  With  the  knowledge  of  the  channel  coefficients  at  the  ith 
node,  a  matched  filtering  operation  is  applied  on  each  of  the  received 
signals  in  the  form  of  ^PSjhj  i/No)  yj Therefore,  the  SNR 
at  the  output  of  the  matched-filter  is  expressed  as 


7 j,i  — 


Psj  \  hj,i 


Nn 


(5) 


After  each  source  Si  has  decoded  its  (TV  —  1)  received  symbols 
{yj,i}f=ij^i’  the  set  of  available  decoded  data  symbols  at  each  node 
in  the  network  is  given  by  the  matrix 
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where  2^  acts  as  a  binary  indicator  function  as  follows 


T  •  _ 

Am 


if  node  Si  decodes  Xj  correctly 
otherwise 


Si 

s2 


Sat-i 

Sat 
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B.  Cooperation  Phase 

The  analysis  of  the  Cooperation  Phase  is  considered  in  this  subsec¬ 
tion,  with  the  assumption  that  the  nodes  are  perfectly  synchronized 
by  a  distributed  algorithm.  In  the  Cooperation  Phase,  each  node  Si 
in  its  assigned  time-slot  receives  a  signal  of  the  N  —  1  source  nodes. 
In  particular,  each  node  Sj  other  than  the  destination  node  S%  forms 
a  linearly-coded  signal  of  the  other  source  nodes’  received  symbols 
and  transmits  it  to  node  Si  during  the  ith  time-slot,  simultaneously. 
Moreover,  at  each  node  Sj,  each  linearly-coded  signal  contains  at 
most  the  received  data  symbols  from  (TV  —  2),  since  each  node 
aims  at  relaying  the  remaining  (TV  —  2)  nodes’  symbols  to  node 
Si.  Specifically,  the  signal  transmitted  to  node  Si  from  node  Sj  is 
composed  from  the  received  data  symbols  of  the  jth  row  in  the 
matrix  X  excluding  the  data  symbol  XiXij  (since  that  symbol  xi 
was  originally  generated  at  node  Si)  as  follows 

Xj  {XkPk,j}k—l3k^j  \  (8) 


Based  on  equation  (6),  during  the  ith  time-slot,  the  signal  transmitted 
from  the  jth  node  is  given  by 
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where  Sk{t)  is  the  signature  waveform  of  the  kth  node  and  Pkj  is 
the  power  at  jth  node  used  to  transmit  the  symbol  Xk-  Moreover,  hj ^ 
is  the  channel  coefficient  between  source  nodes  Sj  and  Sp,  which 
has  already  been  estimated  in  the  Broadcasting  Phase  during  the 
jth  time-slot.  Clearly,  the  functions  f(Xj)  at  each  node  are  linear 
combinations  of  symbols  received  from  other  nodes.  The  operation 
of  the  Cooperation  Phase  can  be  expressed  in  matrix  form  as  follows 


Si 
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Sn 

T\ 

'  0 

■  f(X?) 

X 
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0 

*  *  f{Xh) 
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-f(XiN)  ■ 

■■  mN)  ■ 
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The  received  signal  at  the  ith  node  during  the  ith  time-slot  from 
the  (TV  —  1)  other  nodes  is  given  by 

N  N 

^  ^  ^rn,if{,Xrn)(t)  T  Wj(t) —  ^  ^  T 

rri—1  m= 1 

m^i  m^i 

(11) 

where  Wi(t)  is  the  additive  white  Gaussian  noise  at  node  Si  and 

N 

^i,m  —  E  | hk,i  |  \/ Pfn,  k-^m,ky  (12) 

k= 1 

k^i,k^m 

where  the  summation  in  (12)  contains  at  most  (TV  —  2)  terms, 
depending  on  how  many  data  symbols  have  been  decoded  correctly. 

III.  Multi-source  Signal  Detection 

Based  on  the  received  signal  3^(£),  node  Si  performs  a  multi¬ 
source  detection  operation  to  extract  the  (TV  —  1)  symbols  of  the 
other  nodes.  Each  soft  symbol  Xj,  j  £  {1,2,...,  TV}^-  is  detected 
by  passing  the  received  signal  through  a  Matched  Filter  Bank 
(MFB)  of  (TV  —  1)  branches,  matched  to  the  corresponding  set  of 
the  nodes’  signature  waveforms  sj(t)  for  j  £  {1,2,...,  TV}^-  and 
sampling  them  at  the  end  of  the  symbol  duration  Ts  to  obtain 

N 

^i,j  =  (Ti(^))  sj  (0)  =  ^  ^  xmai,mPm,j  T  wi,ji  (13) 

m= 1 
m^i 


where  pmj  is  the  correlation  coefficient  between  the  signature 
waveforms  sm(t)  and  Sj(t).  The  matched- filtered  signal  forms  an 
(TV  —  1)  x  1  vector  comprising  all  the  3^,/s  signals  as 


y%  —  RiAiXi  +  Wi , 

(14) 

where 

yt  =  [yi,u  ■  ■  ■ ,  , 

(15) 
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(16) 

Wi  =  [w»,  1, . . . ,  w*,i- 1,  •  •  ■ ,  Wi,N\  ~  CAf(O,N0Ri),  with 


1 

Pl,(i-1) 

pi,(i+i) 

Pl,N 

Ri  = 

P(i-  1),1 

1 

P(i-i),(*+i)  • 

1 

’  P(i-l),IV 

P(i+l),l 

•  •  •  P(i+l),(i-l) 

'  ‘  P(i+l),IV 

Pn,i 

PN,(i~  1) 

PN,(i+l) 

1 

(17) 

and 


"fli,  1 

0 

0 

0 

0 

1) 

0 

0 

0 

0 

1) 

0 

_  0 

0 

0 

where  both  matrices  Ri  and  Ai  have  dimensions  (TV  —  1)  x  (TV  —  1). 
The  signal  vector  3^  can  then  be  decorrelated  with  the  assumption 
that  matrix  Ri  is  invertible  with  the  inverse  matrix  R u1  which  yields 


yz  =  Ri1yz=  A^i+Wi,  (i9) 

where  ti,  ~  CAf(0,  A^o-R”1).  Therefore,  the  soft  symbol  of  x3  can 
be  obtained  from  detected  signal  vector  3^  at  node  Si  at  the  output 
of  the  jth  branch  of  the  MFB  and  is  expressed  as 

y%j  dijXj  +  Wi,j,  (20) 


where  Wij  ~  CN{ 0,  TV0ttj)  with  r^j  being  the  jth  diagonal  element 
of  the  inverse  matrix  R~^  associated  with  the  data  symbol  Xj  and 
dij  is  given  by  (12).  Without  loss  of  generality,  let  pj^i  =  p  for  all 
j  /  i.  Thus,  it  can  be  easily  shown  that 

1  +  (TV  —  3)p 


1  +  (TV  —  3)p  —  (TV  —  2  )p2 


:  r N-i. 


(21) 


It  should  be  noted  that  upon  the  completion  of  the  Broadcasting 
Cooperation  Phases,  the  jth  data  symbol  Xj,  j  £  {1,2,...,  TV}^-  is 
relayed  at  most  (TV  —  2)  times  before  reaching  node  Si.  In  addition, 
since  in  the  Broadcasting  Phase,  the  source  node  Sj  has  already 
broadcasted  its  data  symbol  Xj  to  all  the  other  nodes  including  Si, 
this  implies  that  node  Si  received  at  most  a  total  of  (TV  —  1)  copies 
of  the  symbol  xj.  This  fact  will  be  used  later  to  prove  that  a  full 
diversity  order  of  (TV  —  1)  per  node  is  achieved  at  a  high  SNR. 

On  the  other  hand,  since  Tjj  £  {0, 1}  for  i  ^  j,  represents  the 
detection  state  at  Si  of  the  data  symbol  xj\  then,  in  general,  node 
Sk  for  k  £  {1,2,...,  N}k^i,k^j,  forwards  the  symbol  xj  to  node 
Si  only  if  it  has  successfully  detected  it.  Therefore,  at  node  Si,  all 
the  Tj,k$  form  a  binary  (base-2)  number 

ij  —  •  •  •  Rj,k=(i- 1-1)  3'j,k=(i—l)  •  •  •  -22/ ,  1  ]  2  \  -Ej,j  7  (22) 


that  represents  one  of  the  2^N~2^  detection  states  of  the  (TV  —  2) 
nodes  Sj9 s  acting  as  relay  nodes.  In  other  words,  |Z}|  contains  at 
most  (TV  —  2)  l’s.  For  example,  for  TV  =  4,  the  detection  state  of  x\ 
at  node  S2  is  given  by  1 \  —  \T\^  2’1j3]2  with  \l\\  taking  decimal 
values  in  {0, 1,2,3}  in  the  form  of  Z\  —  [0  0)2,  [0  1)2,  [1  0)2,  or 
[1  1]  2  ^  respectively. 

The  detection  of  data  symbol  xj  at  the  node  Si  can  be  achieved 
through  combining  the  signals  received  in  the  Broadcasting  and 
Cooperation  Phases.  However,  during  the  Cooperation  Phase,  it 
might  occur  that  aij  =0  (i.e.  the  receiving  node  might  not  receive 
any  linear  combination  for  some  symbol  Xj  and  this  occurs  when  all 
the  other  nodes  decode  Xj  incorrectly,  simultaneously).  This  in  turn 
implies  that  the  detection  state  for  symbol  xj  at  node  Si  is  |Zj|  —  0. 
Thus,  the  jth  branch  of  the  MFB  j)i,j  is  not  added  into  the  detected 
Xj.  Hence,  for  |Zj|  =  0,  a  simple  phase  correction  is  applied  to 
the  received  signal  yjj  during  the  jth  time-slot  in  the  Broadcasting 
Phase  by  multiplying  it  with  the  deterministic  conjugate  of  the  known 
channel  phase  response  (i.e.  Xj  =  e~^^iyj^).  However,  for  |Zj|  >  0 
(i.e.  dij  ^  0),  the  phase-corrected  signal  is  also  combined  with  yi:j 
which  yields 
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x)  = .  +  yid 


(23) 


where  yjj  =  \JPs3 hj^Xj  +  rijp.  Upon  substitution  of  yj^,  (12)  and 
(20)  into  (23),  the  combined  symbol  x3  at  node  Si  becomes 

(  N  \ 


^  ^  I hk,i  |  yj Pj,kPj,k 


+  (24) 


V 


fc=l 


J 


where  rjij  —  e 


—  P  j4>j 


+  Wij  is  the  zero-mean  equivalent  noise 


with  variance  (1  +  rjv-i)TVo,  with  the  assumption  that  n3 and  Wjj 
are  statistically  independent.  Hence,  the  conditional  signal-to-noise 
ratio  given  the  detection  state  Xj  of  the  combined  data  symbol  xj  at 
node  Si  can  be  expressed  as 

I  hj,i\2PSj 


N0  ’ 

=  \  (  '-,.k/'v.>:va.=i  i ^ 

k^i,k^j 


(l+rjv-i)iVo 


if  |2j|=0 


if  \TA  >  0 


(25) 


The  total  transmit  power  Pj  associated  with  transmitting  symbol 
Xj  is  distributed  among  the  (TV  —  1)  transmissions.  In  particular, 
Pj  =  P8  ■  -\-V]N k=i  Pi  k,  where  Ps  is  the  transmit  power  at  node 

J  ^k^i.kfj  3 

Sj ;  while  the  Pj,k  s  are  the  power  allocations  for  the  remaining 
transmissions  at  the  other  nodes.  Without  any  loss  of  generality,  it  is 
assumed  that  all  the  symbols  are  assigned  the  same  total  power,  such 
that  P  =  Pj,\/j  G  {1,2,...,  TV}.  For  simplicity,  it  is  further  assumed 
that  Ps  =  PSj  =  Pj:k  m  P/{N  —  1),  Vfc  G  {1,2, .. . ,  TV}^^^-. 
However,  it  should  be  noted  that  optimal  power  allocation  among 
the  different  nodes  is  also  possible;  however  is  beyond  the  scope  of 
this  paper  as  it  would  deviate  from  the  main  thrust  of  this  work. 


IV.  Symbol  Error  Rate  (SER)  Performance  Analysis 


In  this  section,  the  exact  SER  expression  for  the  M-ary  Phase 
Shift  Keying  (M-PSK)  modulation  for  the  DF  protocol  of  the 
symbol  xj  detected  at  node  Si  is  derived.  Since  the  detection  at 
each  node  is  statistically  independent  from  the  others,  X^’s  for 
k  G  {1,  2, . . . ,  N}k^i^j  are  defined  as  independent  Bernoulli 
random  variables  with  a  distribution  expressed  as  [4] 


T(I ,-fc) 


1  —  Pj,ki  if  Pj.k  —  1 

Pj,k,  if  Pj,k  —  0  5 


(26) 


where  Vj^  is  the  SER  of  detecting  Xj  at  the  Sk  node.  Thus,  the 
probability  of  Xj  detection  in  state  2}  at  node  Si  is  written  as 

N 

Pr(|2}|)=  J]  (27) 

k  =  l 

k^i,k^j 


where  in  the  example  of  TV  =  4,  for  |X2|  —  1  (i.e.  X2  —  [0  1]2), 
Pr  (\I2\  —  l)  =  Pi, 4  (1  —  Pi, 3).  In  general,  the  conditional  SER  of 
M-PSK  systems  with  the  instantaneous  signal-to-noise  ratio  (SNR) 
7  given  a  generic  set  of  channel  coefficients  {/z}  is  expressed  as  [10] 

KM-D./M  (  bpsK^ 

JO  6XP\  sin2  9  ) 


- 1  [ 
77  J  0 


d6,  (28) 


where  bpsx  —  sin2(7r/M).  Based  on  (4  -  5),  the  SNR  in  detecting 
the  symbol  xj  at  node  Sk  given  the  channel  gain  is  7 ^  = 
PSj  \hj^k  |2/TVo.  In  general,  the  magnitude  squared  of  a  circularly  sym¬ 
metric  Gaussian  random  variable  hj ^  is  modeled  as  an  exponential 
random  variable  with  rate  1/cr2 (i.e.  \h3^\2  ~  Sxp(l/a2k)),  where 


a2k  is  the  channel  gain.  Thus,  by  averaging  the  expression  in  (28) 
with  respect  to  \h3^\2,  the  SER  of  detecting  Xj  at  Sk  is  obtained  as 


Pj,k  — 


K  Jo 


(M~'P/M  ^fbPSKPs 


G 


\No  £ 


2 

1  OjX 


dO , 


(29) 


where  G(w(0),  a2)  is  defined  in  (54)  in  Appendix  I.  On  the  other 
hand,  based  on  the  conditional  SNR  7 ^  expression  in  (24),  the 

conditional  SER  of  symbol  xj  at  node  Si  can  be  shown  to  be 


(7'izi)  =  - 


l 


exp 


(  bpSK^j\x)\ 

y  sin2  9  J 


d6 ,  (30) 


which  reduces  to  the  following  two  cases: 

1)  |Xj  |  =  0:  After  averaging  over  the  exponential  random  variable 


\hj,i\2,  the  conditional  SER  can  be  expressed  as 


^(7: 


j  \\P\=oJ 


K  Jo 


(M  —  1)tt/M  /  bpSKps 

\  TVq  sin2  9  ’ 


dJ9.  (31) 


2)  m  >  0:  Similarly,  the  conditional  SER  can  be  obtained  as 


if 

*  Jo 


exp 


PSK  0/2 


where 


and 


2  sin2  9 
2 


XCi,j  —  \hj,i\  T  y  ^  \hk,i\Xj,k 


k=  1 

k^i,k=£j 


CpSK  = 


ZbpsxPs 
(1  +  VN-i)Nq 


d9 , 

(32) 

(33) 


(34) 


Clearly,  H2-  is  a  sum  of  1  +  |Xj|  Rayleigh  random  variables. 
The  analysis  for  the  conditional  SER  is  analogous  to  that  of 
the  equal  gain  combining  in  [11],  using  the  Gauss-Hermite 
formula  [12,  p.  890,  eq. (25 .4.46)].  Thus,  after  averaging  over 
the  channel  statistics,  the  conditional  SER  is  given  by  [11] 

r(M-l)7r/M  ^ 


^(7]||Tj|>o)  -  2?r2 


—  f 

K2  Jo 


7^1 


wniP j 


,  9  Jd9, 
(35) 


where  /cn,  wn  are  the  zeros  and  weight  factors  as  given  in 
[12,  p.  924,  table  (25.10)]  and  Np  is  the  order  of  the  Hermite 
polynomial  Hnp(-),  respectively.  It  was  verified  that  Np  —  20 
results  in  excellent  accuracy.  In  addition  [11] 


iVj(0) 


sin2  (9)  1 
2  Cpsk  4 


/ 


ali  + 


V 


N 

E 

k= 1 

k^i,k^j 


O'?  hXi. 


Jj,kxj,k 


(36) 


and  ifFj{ v,  0)  =  iRj(p,  0)  cos  (iOj( v,  0)),  with 
iRj(  v,0)  =  ^/X*{6)+Y*{v,  61)  x 

D 


(37) 


k= 1 

ky^i,k^j 


and  iSj(u,0)  is  defined  in  (38)  (top  of  next  page).  On  the 
other  hand,  X(0),  Y(is,9)  are  defined  as 


X[0) 


/it  sin (9) 

2  Cpsk  ’ 


and  Y (v,  9) 


vsin 2  (9)  (\  3  sin2(6>)\ 

Cpsk  1  \2’  2’  2 Cpsk  J 


respectively,  where  sgn(- )  is  the  sign  function  and  \Fi  (•;•;•) 
is  the  Kummer  confluent  hypergeometric  function  [12].  More¬ 
over,  A(v,t)  and  B(y,r)  are  defined  as 
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0)  =  arctan 


(Y{y,9)\ 
\  X(9)  ) 


+  arctan 


U  (y,°id) 


N 

^  arctan 
k= 1 

ky£i,k^j 


(  B{v,  ffjkpk)  \ 


+  f  N  ~  sgn  {Y(v,  9))  -  sgn  (B(u,  cr^))  -  ^  sgn  (B(v,  a2^kXjtk)) 


k= 1 

k^i,k^j 


(38) 


\p)  [^PSKa2(1  +  Ke) 


e+i  F1  (^  +  |, \,l  +  l ,£+ 


2  ( 

V  M  I 


(l+^*iV —  1  ^  Nq 
^PbpsK^ii  1+f  f) 


sin2((M^; 


tt{2£  +  3) 


(M  —  1)-7T 

M 


j2£+3 

(50) 


=  iFi  ^-1;  1;  ,  and  B(v,t)  =  T  \M 

respectively;  with  T(-)  being  the  gamma  function  [12].  Also, 
Dj  i  is  defined  as 

Dp  =  (4D 

Given  the  detection  state  Xj,  which  takes  2^N~2^  values,  the  SER 
for  detecting  the  data  symbol  xj  at  the  ith  node  can  be  calculated 
using  the  law  of  total  probability  as 

2(iV-2)_1 

P3ser  =  xj  I  \Zj  I;  =  '  Pr(l^l  —  -0 

£=0 

2(iV  — 2)  _  i 

=  'f(7;nx.|=o)-Pr(|i;i  =  0)+  £  ^(7lMx*.|=f)  '  Pr(|7j|  =  £) 

3  1=1  3 

(42) 

where  Pr(xj  ^  xj\l))  is  the  probability  of  making  a  symbol  error 
for  a  particular  detection  state,  and  Pr(|Xj|)  is  as  defined  (27). 

V.  Asymptotic  Upper  Bound  SER  Analysis 

The  asymptotic  upper-bound  SER  performance  is  obtained  at  a 
high  SNR  by  performing  a  series  of  approximations  to  the  term 
f°r  ^  £  {0,1,...,  2^~2)  —  1}  and  also  the  term 

Pr(|Xj|).  It  should  be  noted  that  finding  the  distribution  of  the 
sum  of  independent  but  not  identical  Rayleigh  random  variables  as 
given  by  7i23  in  (33)  can  be  extremely  difficult.  Therefore,  it  is 
assumed  that  the  channel  gains  are  identical  (i.e.  a2 ^  —  cr2  i  — 
o 2,  Wk  £  {1,  2, . . . ,  N}k^i).  This  implies  that  the  same  average 
SNR/symbol/path  is  assumed  between  the  nodes  in  the  network. 

At  high  SNR,  it  is  expected  that  the  SER  Vj^  of  detecting  xj 
at  node  Sk  for  k  £  {1,  2, . .  * ,  becomes  sufficiently 

small  such  that  1  —  Pj  k  ~  1.  Thus,  only  the  terms  in  the 
quantity  Pr(|X*|)  =  n  &=i  Y(X,-  k)  that  will  count  are  those 

corresponding  to  the  nodes  that  have  decoded  their  received  Xj 
symbol  incorrectly  [4].  Hence,  let  o  and  i  <3>j  denote  the  sub¬ 
sets  of  the  indices  of  the  nodes  that  decode  Xj  erroneously  and 
correctly,  respectively.  That  is,  o^j  =  {k  •  Pj,k  =  0}  and 
l =  {k  :  lj,k  =  1},  for  k  £  {1,  2, . . . ,  N}k^k^j.  Moreover, 
\o®j\  and  |i^j|  £  {0, 1, . . . ,  N  —  2}.  In  addition,  it  should  be 
noted  that  |0^j|  +  |i^j|  =  (N  —  2)  for  any  detection  state  Xj. 
Also,  since,  the  transmitter  powers  allocated  so  source  node  Sj  and 
the  remaining  nodes  for  the  transmission  of  symbol  Xj  is  defined 
as  a  ratio  of  the  total  power  P  as  Ps  —  PSj  —  Pj ^  =  /iP, 
\/k  £  {1,  2, . . . ,  N}k^i,k^j  with  ii  —  l/(N  —  1),  the  expression 


for  Pr(|Xj|)  given  in  (27)  can  then  be  expressed  with  the  aid  of  (55) 
(see  Appendix  I)  as 

N  /  AT  \  \0$j\  /  /rt  \  |0^j| 

,  (43) 


Pr(Pj|)  —  n  T(A4f 

k  (E  o  j 

where 

^  r(M-l)7r/M 


V  =  ~  , 
K  J  0 


f 


sin 2  (6)  dO  = 


libpsKcr2  y 

27r(M-l)  +  Msin(^) 
47 tM 


(44) 

With  respect  to  ^(7j||xq=o)’  uPPer_bound  approximation  at  high 
SNR  can  be  expressed  with  the  aid  of  (55)  as 


^(7*||ii|=o)  -  F 


(  bpsxPs 


< 


N0 

P 


flbpSK°i 


(45) 


V  A"o  sin2  0 1 

In  order  to  find  an  upper-bound  approximation  for  ^(yj^i^^),  the 
expression  (33)  must  be  approximated.  Based  on  (25),  the  conditional 
SNR  for  1 1)  |  =  t  is  given  by 

/  \2 

1^1+  X  IMA/.  •  (46) 


AlUjM  0  +  rN-i)No 


V 


k= 1 

ky^i,k^j 


) 


Thus,  the  average  SNR  conditioned  on  the  network  detection  state 
can  be  shown  to  be  [11] 

7av,(M-i)  —  7 si  5  (47) 

where  /ySi  =  Psg2 /(I  +  vn-i)No.  An  accurate  approximation  to  the 
conditional  SER  has  been  determined  in  [11]  and  is  expressed  as 

=  U+l)  (bpSKlav,{i+\)i  - — JpNj  ■  (48) 


In  general 


Therefore,  with  the  aid  of  (56)  (see  Appendix  II)  can 

be  shown  to  be  as  expressed  in  (50)  (top  of  page).  Therefore,  the 
asymptotic  SER  expression,  after  a  series  of  manipulations,  can  be 
written  as  (51)  (top  of  next  page),  where  it  should  be  noted  that 

1  •  ki-,i-i,<i>,i7(-v-  D- 

The  cooperative  diversity  order  of  a  wireless  system  is  identified 
from  the  SER  expression  as  follows  [4] 


Her  ~  ( SNR  ■  A)~d ,  (52) 

where  SNR  =  P/No  is  the  signal-to-noise  ratio  term,  the  exponent  d 
denotes  the  diversity  order  and  A  defines  the  cooperation  gain.  Thus, 
it  is  clear  that  d  =  N  —  1  and  the  STNC  achieves  full  diversity. 
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ll*il+1  2F1  (li^l  +  j  |i*j|  +  |,sm2 

7r(^PSi^crf)iV_1(2|l^>jl  +  3) 
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Fig.  2.  QPSK  SER  Performance  of  the  STNC  for  N  =  3  and  N  =  4  Nodes 


VI.  SER  Performance  Evaluation 


G(w(d),o*)  =  Jo  exp  (-«(%)  (E)  exp  (-£)  dy  =  ±  +  J(g)g2  , 

(54) 

where  w{6)  is  some  function  of  0. 

For  large  values  of  w(0),  the  denominator  of  G(w(0),  a2)  can  be 
approximated  as  1  +  w(0)a2  «  w(0)a2.  Thus, 

G{w{0),a2)  <  ^^2-  (55) 

IX.  Appendix  II 

The  term  Ip  (£,  (j>)  can  be  shown  to  be 

L 


II 


1  G 

«>*)  =  £  / 

^  Jo 


K  sin2  f>  +  £ 

(l/£)L(sin(<«)2L+1 


df> 


sin2(0) 


(56) 


where 


In  this  section,  the  simulated  QPSK  SER  performance  and  the 
theoretical  exact  and  upper-bound  SER  expressions  of  the  STNC 
scheme  with  N  —  3  and  N  —  4  nodes  are  evaluated.  In  the  following 
simulations,  all  the  channel  coefficients  are  i.i.d.  with  equal  unity 
gain  (i.e.  hi:j  +-  CM(0, 1),  Vi,  j  €  {1, 2, . . . ,  N}  and  i  j). 

It  is  clear  from  Fig.  2,  that  as  the  number  of  nodes  increases, 
the  better  the  performance.  This  is  due  to  the  fact  that  with  the 
increase  in  N,  higher  diversity  gains  are  achieved.  In  addition,  it 
is  noticed  that  the  best  SER  performance  is  achieved  when  p  —  0, 
that  is  the  signature  waveforms  are  perfectly  orthogonal  and  there 
is  no  interference  between  the  nodes.  However,  when  p  —  0.75, 
there  is  about  2  dB  SER  performance  degradation.  This  in  turn 
signifies  the  importance  of  the  STNC  scheme  as  it  allows  N  nodes  to 
communicate  simultaneously  and  achieves  N  —  1  diversity  order  with 
only  a  slight  degradation  for  non-orthogonal  signature  waveforms. 

VII.  Conclusions 

In  this  paper,  the  novel  Space-Time  Network  Coding  (STNC) 
scheme  that  allows  N  nodes  to  exchange  their  data  symbols  over 
a  total  of  2N  time- slots,  was  presented.  The  exact  and  asymptotic 
SER  expressions  were  derived  and  it  was  shown  that  with  the  STNC 
scheme,  each  node  achieves  a  full  diversity  order  of  N  —  1.  Finally, 
it  is  concluded  that  the  STNC  serves  as  a  potential  many-to-many 
cooperative  communication  scheme  and  its  scope  goes  much  further 
beyond  the  generic  source-relay-destination  communications. 

VIII.  Appendix  I 

Let  Y  be  an  exponentially  distributed  random  variable  with 
probability  density  function 

fv(y )  =  ~2  exP  (“  J2)  >  y  -  (53) 

where  a2  is  the  rate  parameter.  Averaging  the  function  exp(— w(6)y) 
over  the  distribution  of  the  exponential  random  variable  Y,  yields 


*i(a;/3,/3';C;*,2/>=  X  X 

m— 0 n= 0 


{a)m+n((3)m((3  )n  m  n 

ra!n!(C)m+n 


(57) 


is  the  Appell  hypergeometric  function  of  the  first  kind  [13].  It  can 
be  shown  that  for  y  — ►  0,  the  Appell  hypergeometric  function  can 
be  approximated  as  [13] 


F1  (a;  /?,  ff ;  C;  x,  y)  «  2F\  O,  P,  75  x) ,  (58) 


where  2Fi  (<t5/^;7;^)  is  the  hypergeometric  function  [12]. 
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